Atomistic theory of the shear band direction in amorphous sohds 
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One of the major theoretical riddles in shear banding instabilities is the angle that the shear band 
chooses spontaneously with respect to the principal stress axis. Here we employ our recent atomistic 
theory to compute analytically the angle in terms of the characteristics of the Eshelby inclusion that 
models faithfully the eigenfunction of the Hessian matrix that goes soft at the plastic instability. We 
show that loading protocols that do not conserve volume result in shear bands at angles different 
from 45° to the strain axis; only when the external strains preserve volume like in pure shear, the 
shear bands align precisely at 45° to the strain axis. We compute an analytic formula for the angle 
of the shear band in terms of the characteristics of the loading protocol; quantitative agreement 
with computer simulations is demonstrated. 
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FIG. 1. Color Online: The shear band that occurs in a 2- 
dimensional amorphous solid upon uniaxial compression (up- 
per panel) and extension (lower panel). The angle with re- 
spect to the principal stress component measured in the upper 
panel is 46° ± l'', whereas in the lower panel it is 54° ± l''. 
This Letter provides a microscopic theory of this asymmetry. 

The riddle that this Letter aims to solve is demon- 
strated in Fig. [TJ Amorphous solids put under different 
loading conditions display shear bands that appear at 
different angles with respect to the principal stress axis 
[l|. In past workphenomenological models were applied 
to this riddle [2|-|j] , but a microscopic approach was lack- 



ing. In this Letter we provide a microscopic theory of 
this phenomenon in athermal systems that are strained 
under quasi-static conditions [10. The Letter culmi- 
nates with an analytic formula for the angle in terms of 
one parameter which is determined by the loading condi- 
tions. For simplicity and concreteness we will explain the 
theory in 2-diniensions, the generalization to 3-dimension 
is available [7|. 

Recently there had been rapid progress in under- 
standing the nature of plastic instabilities in Ather- 
mal Quasistatic (AQS) conditions |8l4l2|. Consider a 
glass-forming system made of TV particles interacting via 
generic potentials (i.e. at least twice differcntiable ev- 
erywhere). The total energy can be written at T = 
in terms of the positions ri, r2, • • • tat of these particles, 
U = U{ri,r2- ■ • , Tat). The Hessian matrix is defined as 
the second derivative llSl 



Hi 



a^[/(ri,r2--- ,rjv) 
dridrj 



(1) 



The Hessian is real and symmetric, and therefore can be 
diagonalized. Excluding Goldstone modes whose eigen- 
values are zero due to continuous symmetries, all the 
other eigenvalues are real and positive as long as the sys- 
tem is mechanically stable. In equilibrium, without any 
mechanical loading, eigenfunctions associated with large 
eigenvalues are localized due to Anderson localization. 
But all the eigenfunctions associated with low eigenval- 
ues (including all the excess modes that are typical to 
amorphous solids) are extended. At low values of exter- 
nal loading one observes "fundamental plastic instabili- 
ties" when eigenvalues of some modes approach zero via 



a saddle node bifurcation [IJ]. Simultaneously the as- 



sociated eigenfunction localizes on a typical quadrupolar 
structure (in both 2 and 3 dimensions) that is identi- 
cal to the non-afSne displacement field associated with 
the plastic event. An example of such a localized eigen- 
function that is seen upon uniaxial compression is shown 
in Fig. [21 At higher values of the external strain the 
nature of the plastic instability can change drastically. 
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FIG. 2. Upper panel: A fundamental plastic event during 
an AQS uniaxial compression simulation of a 2-dimensional 
amorphous solid. Shown is the non-afiine displacement field 
in the whole system. Lower panel: a small window around 
the core of the event shown in the upper panel. For clarity 
we show only the incoming and outgoing arrows. 

I 



For well quenched amorphous solids one finds that in- 
stead of single fundamental events there exists a bifurca- 
tion to the simultaneous appearance of a finite density of 
such events, (i.e. infinity of them in the thermodynamic 
limit), and that in tvifo dimensions they organize along a 



line [llj, |12[ . In previous work the angle of this line with 
respect to the principal stress axis was calculated analyt- 
ically for external shear, with the result of 45°. We will 
explain now why external shear is special, and why the 
angle will change for other loading conditions. 

A very important realization is that the fundamental 
plastic event can be excellently modeled analytically as 
an Eshelby inclusion [15| . As is well known, the Eshelby 
construction in 2-dimensions starts with a circle taken 
out from an elastic medium, made into an ellipse that is 
then pushed back to the medium. This is the Eshelby 
inclusion. In terms of the two orthogonal directions n 
and k one writes the eigenstrain of the inclusion as 



ta^ 



Cnflanfj + Ckkakp . 



(2) 



The "eigenvalues" Cn and Cfc characterize the relative 
magnitude of the non-affine displacement associated with 
the fundamental plastic instability. The direction h is 
that of the principal stress axis, and k the orthogonal one. 
The ratio |Cra/Cfc| depends, as shown below, on the load- 
ing conditions. When area preserving prevails, Cn = ~Cfe 
111 |l2|, but in general these numbers are independent, 
and their ratio characterizes the loading condition. We 
argue below that this ratio determines uniquely the angle 
of the shear band once there appears a density of such 
events with a prescribed ratio |Cn/Cfe|- 

Having the eigenstrain ([2]) one asks what is the dis- 
placement field associated with such an inclusion. The 
calculation is shown in detail in the accompanying mate- 
rial [l™, with the final result 
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with X denoting an arbitrary Cartesian point in the ma- 
terial, r = \X\, f = X/r and a is a parameter known as 
the "core size" . The stress field induced by this displace- 
ment field is denoted below as cr^a- 

As said above under simple shear one conserves area. 



and C„ = — Cfe. Thus the first term inside the square 
parenthesis in Eq. [3] vanishes, and in the two other terms 
Cn — Cfc = 2Cn = e, cf. [12| . We have seen in Ref. [12] , and 
we will see below, that when such a condition applies, the 
angle of the shear band is precisely 45° with respect to 



the principal stress axis. But for the more general loading 
conditions in which C„ ^ —C,k the angle is different, see 
Eq. dH]) below. 

When the amorphous solid is well quenched, the nature 
of the plastic instability changes when the external strain 
becomes sufficiently high; a density of inclusions of the 
form ([3]) appear simultaneously. To find their geometrical 
organization we need to compute the energy of such M 
inclusion in a system of total volume V and minimize it 
with respect to their orientation and geometry. The total 
energy of such M inclusion has four contributions: 



E — Ejnat + -E-cxj + Eesh + Einc 

with each component of energy defined as 
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(5) 
(6) 

(7) 

(8) 
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Here the eigen-strain e^l and volume Vq associated with 
any i*'* Eshelby inclusion are given as 



yj = Tra" 



^afi 



o [2n^np - Sap) + da^ (9) 



The reader should recognize that the first contribution, 
^, is the energy of the elastic matrix and the second 
(([HI)) is the energy due to the interaction of the elastic 
matrix with the inclusions. The third, ^, is the self- 
energy of all the inclusions, and the last, ([8|), is the in- 
teraction energy of different inclusions. 

Also for a 2D material being loaded under uni-axial 
strain with free boundaries along y, we can write the 
form of the global stress tensor as 







(10) 



By Hookc's law, we get the expression for applied global 
stress tensor using the Lame coefficients A and fi: 



Taking trace of Eq. ()lip , we find 

„ 1 ^ 1 



Plugging Eq. ^ in Eq. ^, we find 

00 _ 4a^(A + A^)7 
X + 2^i 



(11) 



(12) 



(13) 



where 7 is the external strain. 

Clearly, the first contribution Eq. ([5]) is independent of 
the distribution of inclusions and will not affect the state 
of minimal energy. The second term Eq. ^ is the only 
one proportional to the external strain 7, and therefore 
for sufficiently large 7 it needs to be minimized first. Be- 
sides 7, this term depends only on the orientation of each 
inclusion, i.e. on the angle </> = cos~"'^(n^) = sin~^(n' ). 
Minimizing the term with respect to this angle we find 
that the minimum is obtained when = or 7r/2. In 
other words, each inclusion is oriented with one axis par- 
allel and the other orthogonal to the uniaxial direction. 
Plugging this information into Eqs. ([7]) and ([5|) simplifies 
them considerably, and see the supplementary material 
for full details [16|. The organization and orientation 
of the density of inclusions is determined by minimizing 
i^inc of Eq. ([8|) . The minimum energy is obtained when 
all the inclusions align on a line that is at an angle 9 with 
respect to the uniaxial direction [l6| . The final prediction 
for this angle is 



cos^ei 



1 



Cn + Cfc 

4(Cn - a) ' 



or 




(14) 



(15) 



One learns immediately from Eq. ([TC)) that the area 
preserving case, C„ = —Cfc, leads to an angle of precisely 
45°. Any other loading condition would result in a differ- 
ent angle. The two extreme cases occur for |Cn/Cfe| ~^ 
and |Cn/Cfc| ^ 00 ■ The first case would realize an angle 
of 30° and the second an angle of 60°. All experimental 
realizations should fall between these two extreme uni- 
versal limits without exception. Indeed, examining the 
wide range of angles reported in Ref. [l[ we see that all 
the data falls within our theoretical limits. We return 
now to our simulations shown in Fig. [l]to rationalize the 
angles observed. 

In order to understand the angle seen in a particular 
experiment one needs to figure how the values of C„ and 
Cfc are determined by the loading conditions. For example 
in uniaxial extension the outgoing direction may become 
highly dominant compared to the incoming one, leading 
to a high ratio of |Cn/Cfc|- This is not the case for com- 
pression, as one can see from Fig. [5J We expect such 
asymmetry to appear in any generic potential due to the 
steep rise of the repulsive part vs the moderate attrac- 
tive tail. The ratio of eigenvalues is obtained from the 
lower panel by averaging the length of the incoming and 
outgoing vectors around the core. We find that the ratio 
of the average length provides a fair estimate of |Cn/Cfc|: 
in this case |Cn/Cfc| ~ 1-15- Using this value in Eq. P^ 
we find an angle 9 « 46° in perfect agreement with the 
observed angle in Fig. [1] upper panel. The exercise is 
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FIG. 3. Upper panel: A fundamental plastic event during an 
AQS uniaxial extension simulation of a 2-dimensional amor- 
phous solid. Shown is the non-afhne displacement field in the 
whole system. Lower panel: a small window around the core 
of the event shown in the upper panel. For clarity we show 
only the incoming and outgoing arrows. 



repeated for uniaxial extension, see Fig. [31 Obviously, in 
the extension case the outgoing arrows are much longer 
than the incoming ones. Averaging the length of the out- 
going and incoming vectors and taking the ratio we find 
for the present case |Cn/Cfc| ~ 4.05. Plugging this value 
into Eq. ((TS|) we get in this case 9 « 54°, again in perfect 
accord with the observed angle in Fig. [T] lower panel. 

To summarize, we note that the atomistic theory is 
based in the understanding that shear localization occurs 
due to a plastic instability in which a density of Eshelby- 
like quadrupolar structures organize on a line. Each in- 
dividual quadrupolar structure is sensitive to the loading 
conditions, parameterized in the 2D case by the ratio 
of "eigenvalues" Cn/Cfc- At zero temperature and under 



quasi-static conditions one is well justified to find the 
preferred geometry by energy minimization. This leads 
to the prediction that all the quadrupoles have the same 
orientation (for sufficiently high external strain) and that 
they organize on a line at an angle 9 with respect to the 
principal stress axis. The angle 9 is limited in this theory 
to fall between 30° and 60°, with 45° being special to the 
area preserving situation Cn/Cfc = ~1- This range of al- 
lowed angles appears to be well in accord with published 
experimental results. 

We have considered in detail the orientation of the 
shear bands under uniaxial compression and extension 
and demonstrated the sensitivity of the ratio of 'eigen- 
values' to these different loading conditions. Finally, we 
showed that the theoretical prediction Eq. (fT5)) is in ex- 
cellent agreement with the observed angles in these load- 
ing conditions. 

The extension of the theory to finite temperatures and 
strain rates is beyond the scope of this Letter, but see 
Ref. [l7[ for a possible direction to go. Also, in the 



follow up paper to this Letter we will present the theory 
for the asymmetry in the yield-strain as a function of C„ 
and Cfc. Stay tuned. 
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